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MONTE-CARLO  SIMULATION  0E-POLARIZATION- SELECTIVE 
SPECTRAL  HOLE  BURNING  IN  FRACTAL  CLUSTERS 

TAPIO  T.  RANTALA MARK  I.  STOCKMAN and  THOMAS  F.  GEORGE ’ 

■  Department  of  Physics.  University  of  Ouiu.  SF-90570.  Ouiu.  Finland 
:  Departments  of  Physics  Astronomy  and  Chemistry,  Center  for  Electronic  fc  Electro- 
ootic  Materials.  239  Fronczak  Hall.  State  University  of  New  York  at  Buffalo.  Buffalo. 

NY  14260 

:  Institute  of  Automation  Electrometry,  Siberian  Branch  of  the  USSR  Academy  of 
Sciences 

INTRODUCTION 

.1  Fractals,  objects  of  noninteger  dimensionality  embedded  into  usual  (three  dimen¬ 
sional)  space,  possess  many  nontrivial  mathematical  properties»4lAi+:  There  exists  a  ^ 
senes  of  physical  systems  whose  geometry  can  be  adequately  described  as  fractal,13-5b 
We  shall  restrict  ourselves  to  the  study  of  one  of  such  physical  realizations  of  fractals, 
namely  fractal  clusters^)  (simply  called  beiow  as  fractals).  More  specifically,  we  shall 
consider  theoretically  the  selective  photomodification  of  fractals  induced  by^-powerful 
iaser  radiation.  Experimentally,  this  phenomenon  has  recently  been  observeato)  in  silver 
clusters  in  colloidal  solutions  and  gels. 

The  effect  of  selective  photomodification  observed  in  Ref.  [61  consists  in  burning 
out  the  spectral  hole  in  the  absorption  contour  of  a  fractal  at  the  frequency  close  to 
that  of  the  exciting  laser  radiation.  The  spectral  hole  exists  at  times  on  the  ordeT  of 
months  or  longer  and,  thus,  can  be  considered  as  persistent.  The  primary  hole  is  highly 
dichroic:  it  is  observed  only  in  the  polarization  of  the  probe  light  coinciding  with  that  of 
the  exciting  radiation  and  is  practically  absent  for  the  normal  polarizations.  In  solutions 
this  dichroism  relaxes  with  the  rotational  diffusion  times  of  the  clusters,  and  in  gels  it 
does  not  significantly  change  for  months.  The  experimental  data'll  indicate  that  the 
photomodification  has  a  threshold  in  the  intensity  of  the  laser  radiation.  Below  we  shall 
exploit  the  last  feature  to  develop  the  theory  of  the  photomodification  which  should 
explain  its  selectivity  and  give  further  predictions  to  stimulate  experimental  study.  [ 

THEORY  AND  NUMERICAL  SIMULATION  OF  SPECTRAL  HOLE  AN 
BURNING  IN  FRACTALS  '  \ 

The  fractal  is  modelled  as  the  set  of  N  polarizable  particles,  called  monomers, 

located  at  the  points  r;,  z  =  1 . N.  The  total  size  of  the  fractal  Rc  is  assumed  to 

be  much  less  that  the  incident  light  wavelength  A.  Therefore  the  electric  field  E0)  of 
the  light  wave  can  be  considered  to  be  the  same  at  the  sites  of  all  the  monomers.  The 
external  plus  local  electric  field  (i.e..  the  integral  field  of  all  the  other  monomers)  induce 
on  the  zth  monomer  a  transition  dipole  moment  d<  oscillating  with  the  light  frequency, 
which  obeys  the  well-known  system  of  equations 

.v 

dxa  =  XoE'*'  -  xo  ^  r~>d]3  .  (1) 

;=  1 


where  the  Greek  subscripts  stand  for  tensor  components  |  summation  over  repeated  in¬ 
dices  is  implied),  and  the  Latin  indices  stand  for  ordinal  numbers  of  monomers;  xo 
is  the  dipolar  polarizability  of  an  individual  (isolated)  monomer:  r,.  =  r.;  -  r;:  and 

n“JI  =  r,,  :r,:. 

The  properties  of  the  monomer  enter  the  system  (1)  oniy  via  the  polarizability 
Xo-  Taking  this  into  account,  we  shall  consider  as  the  spectral  variable  not  the  light 


frequency  u.\  but  the  quantity  A"  =  Rex^1-  As  the  characteristic  of  the  dissipation  in  the 
monomer,  we  shall  use  the  parameter  6  =  Imx^"1 .  To  make  the  connection  to  experiment, 
let  us  specify  the  relation  between  A'  and  u  in  a  model  which  is  realistic,  in  particular, 
for  clusters  obtained  bv  aggregation  in  colloidal  solutions.  This  model  considers  the 
monomer  as  a  macroscopic  sphere  with  radius  Rm ,  consisting  of  the  material  with  a 
complex  dielectric  permittivity  t  =  •'  —  it" .  The  polarizability  of  such  a  sphere  is  given 
by  the  well-known  expression  Xo  =  Rm(e  —  1)U  —  2)-1  .  For  the  metallic  sphere,  the 
permittivity  is  well  described  by  the  Drude  formula  £  =  £0  —  u.-2/iu/(u>  —  17)'  .  where 

r0=const  is  the  contribution  of  the  interband  transitions  (e.g.,  for  silver  to  ~  5),  top  is 
the  electron  plasma  frequency,  and  7  is  the  electron  collision  frequency.  For  most  metals, 
in  particular,  for  the  noble  metals,  7  uiv  and.  consequently,  t"  <C  e' .  In  this  case,  from 
(7)  and  (25),  one  can  see  that  X{w)  turns  to  zero  at  the  point  of  the  surface  plasmon 
resonance,  u>  =  u>3,  where  t1  =  -2.  From  (26)  it  follows  that  ui,  ^  u;p(eo  -r  2)-1/2.  In 
the  vicinity  of  the  surface  plasmon  resonance,  X  and  6  are  simply  expressed  in  terms  of 
the  detuning  from  the  resonance  Q  =  u>  —  u>a  and  the  transition  dipole  moment  of  the 
resonance  1  <£1 2 J  =  3-R^,Au>p/  [2(fo  ~  2)3/,2j ,  namely  X  =  kCl/ldm  -  <5  =  ^7 /2 1 cfi 2 1 

The  theory  is  based  on  the  above  mentioned  assumption  of  the  threshold  character 
of  the  photomodification.  In  addition  it  is  assumed  that  the  photomodification  is  local, 
i.e.  only  those  monomers  are  modified  for  which  the  amplitude  of  the  local  field  exceeds 
the  threshold  magnitude  Eth-  The  latter  depends  on  the  composition  of  the  monomer  and 
the  embedding  medium,  and  also  on  the  specific  mechanism  of  the  photomodification. 
Note  that  the  local  field  is  simply  expressed  through  the  solution  of  Eq.  (1),  E*  = 
dt/xo  •  For  the  sake  of  definiteness,  we  also  suppose  that  the  modified  monomers  do  not 
contribute  to  the  optical  absorption.  This  assumption  is  suggested  by  the  estimates  [6] 
that  the  modification  mechanism,  apparently,  consists  of  melting  and  evaporation  of  the 
monomer  material. 

The  numerical  study  has  been  carried  out  in  the  following  way.  Three  types  of 
fractals  are  generated  and  used  in  the  simulation:  random  walks  (the  fractal  dimension 
D  =  2),  self-avoiding  random  walks  ( D  ~  1.7)  and  Witten-Sander  clusters  ( D  ~  2.5), 
where  the  number  of  clusters  of  each  type  is  1000,  with  mean  number  of  the  monomers 
in  the  cluster  as  N  =  50.  The  results  obtained  for  all  the  three  fractal  types  are  quite 
similar,  so  that  we  shall  give  all  the  numerical  illustrations  only  for  the  random-walk 
fractals. 

For  each  individual  fractal,  the  basic  system  of  equations  (1)  is  solved  numerically 
and  the  factor  Gi  of  enhancement  of  the  local  field  for  the  ith  monomer  is  found,  Gi  — 
|E,|2  /  (E(0))2.  The  above-discussed  condition  of  modification  is  taken  in  the  form 

Gi  >  Gth  ,  Gth=  (EW£<0>)2  .  (2) 

Those  monomers  for  which  this  condition  is  met  are  excluded  from  the  fractal.  After 
that,  the  system  (1)  is  solved  once  again  for  the  modified  fractals  (with  the  excluded 
monomers).  The  optical  absorption  A  per  one  monomer  is  determined  from  the  formula 
A  =  ImXi2,  where  Xzi  =  (dz)  / E^0) ,  and  the  polarization  of  the  exciting  wave  is,  for  the 
sake  of  definiteness,  directed  along  the  z-axis.  Note  that  the  absorption  cross-section  <ra 
per  one  monomer  is  related  to  A  by  aa  =  AnAtZ.  It  can  be  shown  that  the  number  AA 
of  the  modified  (excluded)  monomers  is  expressed  as  AN  =  AA(X)dAf,  where  AA 
is  the  change  of  absorption  upon  the  photomodification. 

In  Ref.  f7l  we  have  found  the  estimate  G  of  the  mean  of  G,  over  the  ensemble  of 
fractals  in  the  scaling  region  (6  <  LYi  <  J?^3)  as 

G  =  (Gi)  ~  Q(R^|AT|)<io+1  . 


(3) 


where  Ro  is  the  characteristic  separation  between  the  nearest  monomers  in  tne  fractal. 

~  0.4  is  the  optical  spectral  dimension,  and  0  =  \  is  the  quality  factor  of  the 

resonance  in  the  monomer.  For  many  monomer  materials,  including  the  noble  metals, 
the  Q-factor  is  large:  e.g..  for  silver.  0  =  3-30  depending  on  the  defect  concentration 
in  the  metal  determined  by  the  preparation  method.  Thus,  the  factor  G  is  large,  which 
means  that  fluctuations  of  the  local  fieid  are  much  grater  than  the  exciting  field  E{0) ■ 
This  fact,  in  the  light  of  Eq.  (2),  allows  us  to  predict  that  the  photomodification  takes 
piace  at  high  values  of  Gf/,  G.  In  the  other  words,  the  characteristic  light  intensity 
of  photomodification  of  fractals  is  much  smaller  (  by  the  factor  ~  G~l  <  1)  than  the 
threshold  intensity  of  the  isolated  monomer  modification.  This  is  a  consequence  of  large 
fluctuations  of  local  fields  in  the  nontrivial  fractals  and.  also,  of  the  high  Q-faccor  of  the 
monomer  resonance  (cf.  Eqs.  (2)  and  (3)). 

In  Fig.  1  we  present  an  example  of  the  dependence  of  the  relative  number  of  modified 
(removed)  monomers  AN/N  on  the  threshold  parameter  Gtk,  with  the  value  X '  of  the 
A'-parameter  for  the  modifying  radiation  as  X'  =  3.0.  One  can  see  that,  in  accord 
with  the  above  arguments,  the  photomodification  starts  at  large  values  of  Gth<  which 
correspond  to  small  exciting  field  intensity  E{0)2 . 
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Fig.  1.  The  relative  amount  of  the 
modified  monomers  as  the  function 
of  the  threshold  parameter  Gth  (2) 
calculated  for  Q  =  3. 


An  example  of  the  calculated  fractal  absorption  contour  is  shown  in  Fig.  2.  One 
can  see  that  this  contour  is  a  broad  peak  centered  at  the  zero  Ar,  corresponding  to 
the  frequencies  near  the  surface-plasmon  resonance.  The  differential  spectrum  of  the 
photoburning,  i.e.,  the  difference  between  the  fractal  absorption  before  and  after  the 
modification,  is  shown  in  Fig.  3  for  the  two  polarizations  of  the  probing  field  and  Xm  = 
-1.  It  follows  from  this  figure  that  for  the  parallel  polarization  the  spectral  hole  is  burned 
out.  centered  at  the  frequency  of  the  exciting  radiation  (X  %  A'*).  For  the  perpendicular 
polarization,  the  spectral  hole  at  X  ==  -Y*  practically  disappears,  while  it  emerged  at 
the  "mirror'’  frequency  at  A’  ~  —A’’’, '2.  This  can  be  understood  in  the  framework  of 
the  binary  approximation  developed  in  Ref.  8.  Thus,  in  fact,  the  spectral  hole  is  highly 
dichroic.  As  the  numerical  results  show,  in  accord  with  the  theoretical  estimates  and 
Fies.  1  and  2.  the  width  of  the  hole  burned  fin  the  A’-variable)  is  on  order  of  Q~l . 
while  the  total  width  of  the  fractal  contour  is  on  order  of  unity  independently  of  the 
Q-parameter.  Thus,  the  total  number  of  holes,  which  can  be  burned  out  in  each  of  the 
two  polarizations,  can  be  estimated  as  approximately  2 Q. 


Fig.  2.  Absorption  spectrum  of  the 
dusters  calculated  for  Q  =  3 


Fig.  3.  Difference  in  the  absorption 
spectrum  caused  by  the  photomodifi¬ 
cation  calculated  for  <5  =  3.  Ar*  =  —  1. 
with  the  polarization  parallel  (curve  1) 
and  perpendicular  (curve  2)  to  that  of 
the  modifying  radiation 


CONCLUDING  DISCUSSION 

The  theory  presented  above  describes  a  qualitative  picture  of  the  polarization- 
selective  persistent  spectral- holes  photoburning  in  fractals,  which  has  been  experimen¬ 
tally  observed  (61.  The  estimate  obtained  on  the  basis  of  the  numerical  results  shows  that 
it  is  possible  to  bum  out  in  a  given  spatial  region  of  a  fractal  media  the  number  of  «  2 Q 
frequency-  and  polarization-resolvable  spectral  holes.  This  number  may  be  on  order  of 
10  for  the  realistic  value  of  Q  —  1  —  10.  Thus,  the  fractal  cluster  media  are  promising 
materials  for  high-density  optical  recording  of  information. 
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